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We report on recent results for the spectrum of heavy quarkonia. Using coarse and anisotropic lattices we 
achieved an unprecedented control over statistical and systematic errors for higher excited states such as exotic 
hybrid states. In a parallel study on isotropic lattices we also investigate the effect of two dynamical flavours on 
the spin structure of charmonium and bottomonium for several symmetric lattices. 



1. INTRODUCTION 

Heavy quarkonia play an important role for 
the theoretical understanding of QCD. Their non- 
relativistic character has frequently been used to 
perform efficient lattice simulations and has trig- 
gered many detailed studies of systematic errors 
such as lattice spacing artefacts, relativistic cor- 
rections and quenching effects. As an additional 
advantage there exists a wealth of experimental 
data, which provide an ultimate check on different 
improvement programmes. Moreover, lattice cal- 
culations have also resulted in predictions from 
first principles for heavy QQg-st&tes containing 
an additional gluonic excitation However, 
those attempts were hampered by the rapidly de- 
caying signal-to-noise ratio of such high-energetic 
hybrid states on conventional lattices. 

More recently anisotropic lattices have been 
used to circumvent this problem by giving the 
lattice a fine temporal resolution whilst maintain- 
ing a coarse discretisation in the spatial direction 
HHI. In a previous study we already reported 
on first NRQCD results for charmonium and bot- 
tomonium hybrid states from anisotropic lattices 



fa] . In Section 2 we report on further applications 
of those methods and study also other excitations 
in heavy quarkonia more carefully. 

In quenched simulations without dynamical sea 
quarks the strong coupling does not run as in na- 
ture and therefore one cannot reproduce experi- 
mental quantities at all scales. Observed devia- 
tions of the quenched hadron spectrum from ex- 
periment have been reported previously and an 
improvement has been noticed once dynamical 
quarks are inserted into the gluon background 
M. Here we study unquenching effects for heavy 
quarkonia and report on our results from isotropic 
lattices in Section 3. 

2. EXCITED QUARKONIA 

To study excited states with small statistical 
errors it is mandatory to have a fine resolution 
in the temporal lattice direction, along which we 
measure the multi-exponential decay of meson 
correlators. To this end we employ an anisotropic 
and spatially coarse gluon action: 
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Here (/3, £) are two parameters, which determine 
the gauge coupling and the anisotropy of the lat- 
tice. Action (Q) is Symanzik- improved and in- 
volves plaquette terms, P^v, as well as rectangles, 
Rfj, v . It is designed to be accurate up to 0(a^, a*), 5T 
classically. To reduce the radiative corrections ^ 
we used self-consistent mean-field improvement 
for both spatial and temporal links. With this 
prescription we expect only small deviations of £ 
from its tree-level value a s /a t . 

For the heavy quark propagation in the gluon 
background we used the NRQCD approach on 
anisotropic lattices as described in ||. From 
the quark propagators we construct meson cor- 
relators for bound states with quantum numbers 
5(0, 1) x L(0, 1, 2) and for hybrid states with addi- 
tional gluonic excitation. For example, the spin- 
singlet operators read 

Q ] Q , Q f A 4 Q , Q^AjAkQ and Q^B t Q . (2) 

Within the NRQCD approach it is paramount to 
establish a scaling region at finite lattice spac- 
ing. Our results in Fig. [j] demonstrate that we 
succeeded in finding such a window. As we only 
measure excitation energies relative to the ground 
state it is natural to present our results as the ra- 
tio R x = (X - IS) /(IP - 15), which gives the 
normalized splitting of state X above the IS. 

For the lowest lying hybrid excitations, ccg and 
bbg, our results from leading order NRQCD |(| 
are in excellent agreement with previous calcu- 
lations on isotropic lattices ^|J3|, but with much 
smaller errors. This is the combined success of 
anisotropic and coarse lattices with a clear sig- 
nal over many timeslices at small computational 
cost. Here we have also checked the spin-averaged 
hybrid against possible finite volume errors, tem- 
poral lattice spacing artefacts and relativistic cor- 
rections, but we did not find any significant effect. 

In Fig. [I] we have also shown our new re- 
sults for higher radial excitations and D-states 
(L = 2). Since all the spin corrections up to 
0(mv 6 ) are now included in our analysis we can 
also determine the spin-structure very accurately. 
In particular we were able to extract the exotic 
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Figure 1. Scaling analysis for excited Quarko- 
nia. We plot the ratio Rx against the spatial 
lattice spacing for different states X = 2S,2P,1D 
and magnetic hybrids (1B,2B). 

hybrid quarkonia, 3 B^ + , explicitly for the first 
time from NRQCD. Our data indicates that the 
spin splittings in hybrid states are enlarged com- 
pared to P-states, whereas the D-state splittings 
are much suppressed. A more detailed discussion 
of the spin structure in quarkonia is presented 
elsewhere ||. 

The dominating systematic error for all the pre- 
dictions in this section is an uncertainty in the 
scale as a result of the quenched approximation. 
This is not yet controlled and we find a varia- 
tion of 10-20%, depending on which experimental 
quantity is used to set the scale. 

3. SPIN STRUCTURE AND SEA- 
QUARK EFFECTS 

To study sea quark effects in heavy quarkonia 
we employ an ensemble of isotropic lattices, where 
the gauge configurations were generated from an 
RG improved gluon action and tadpole-improved 
SW fermions for two flavours of light sea-quarks 
@ . We then performed an NRQCD calculation of 
S- and P-states at two different gauge couplings of 
j3 = 1.80 and — 1.95, corresponding to o« 0.2 
fm and a w 0.15 fm. In addition to the chiral ex- 
trapolation from four different sea quark masses, 
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Figure 2. Chiral extrapolation of the lattice spac- 
ing at two different couplings determined from m p 
(squares) and IP — 15 in Bottomonium (circles). 

we also compared our results directly to quenched 
simulations at the same lattice spacing. 

In Fig. |^ we see an encouraging trend for the 
lattice spacings from different physical quantities 
to agree much better in the chiral limit. This is an 
improvement over quenched simulations, where 
one has large uncertainties in the scale. 

Furthermore, we observed a clear shift upwards 
of the hyperfine splittings as the sea quark mass 
is decreased. For the Bottomonium we find an 
effect of about 3 MeV in a direct comparison to a 
quenched calculation using an identical formula- 
tion of tadpole-improved NRQCD with accuracy 
0(mv 6 ). This is a 10-er effect (or 10%) as shown 
in Fig. ||. In Charmonium the hyperfine splitting 
is also raised by about +15%, i.e. to around 60 
MeV in the chiral limit on our coarsest lattice. 

In P-states one expects a different situation, 
since their wavefunctions vanish at the origin. In- 
deed, the very small hyperfine splitting ( 3 P— 1 Pi) 
is further suppressed (« 3-<r) on dynamical con- 
figurations and there is no resolvable shift for the 
fine structure, e.g. 3 P — 3 Po. This validates the 
quenched approximation for such quantities. 

For the Bottomonium system we are presently 
performing a similar analysis at (3 = 2.1 (a ~ 0.1 
fen) to determine whether those observations still 
hold on finer lattices. 



Figure 3. Spin structure of Bottomonium. Here 
we plot results from quenched and two-flavour 
QCD at the same lattice spacing a « 0.164 fm. 
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